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Abstract. Variability is a defining characteristic of young stellar systems, and optical 
variability has been heavily studied to select and characterize the photospheric proper- 
ties of young stars. In recent years, multi-epoch observations sampling a wider range 
of wavelengths and time-scales have revealed a wealth of time-variable phenomena at 
work during the star formation process. This splinter session was convened to summa- 
rize recent progress in providing improved coverage and understanding of time-variable 
processes in young stars and circumstellar disks. We begin by summarizing results 
from several multi-epoch Spitzer campaigns, which have demonstrated that many young 
stellar objects evidence significant mid-IR variability. While some of these variations 
can be attributed to processes in the stellar photosphere, others appear to trace short 
time-scale changes in the circumstellar disk which can be successfully modeled with 
axisymmetric or non-axisymmetric structures. We also review recent studies probing 
variability at shorter wavelengths that provide evidence for high frequency pulsations 
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associated with accretion outbursts, correlated optical/X-ray variability in Classical T 
Tauri stars, and magnetic reversals in young solar analogs. 



1. Introduction 



The formation and early evolution of stars is inherently a time-domain process: the 
conversion of a dense molecular core into a zero-age main sequence star requires nearly 
every pertinent physical parameter (radius, density, temperature, Vrot, etc.) to change by 
multiple orders of magnitude over the relatively brief timescale of 10s of Myrs. These 
time-scales still dwarf that of a human lifetimq^. and one might a priori conclude that 
star formation is no more amenable to time-domain study than any other aspect of a 
stellar astrophysics. 

Photometric variability has non etheless b een recognized for decades as a com- 
mon trait of many young stars (e.g., |joylll94^ . Historically, this variability has been 
best explored via optical photometry, most often sampling time-scales of days to (a 
few) years. These studies have provi ded de tailed descriptions of the spot properties of 
optically revealed young stars (Vr ba et al.|[T98 8). a comprehensive inventory of stel- 
lar ro tation in young clusters (see first the re view bylHerbst et aLir2007i: Bouvie r et all 
119951 : IStassun et al.iri999l : iRebull et al.l[2004l : ICieza & Baliberll2007l : llrwin et al.ii2008h . 
and a quantitative s tatistical portrait of accretion and extinction-induced variability 
dGrankin et al.ll2007h . These monitoring programs have also identified a num ber of rare. 



astrophysically valuable systems: pre-main sequence eclipsing binaries (e.g., Cargile et al. 



I2OO8. ). disk occulting systems (ie, KH-15D: Hamilton et_al- 2001), and sta rs undergoing 



massive accretion events (i.e., FU Ori, EX Lup, or V1647 - like variables: iHerbig 



1977L 



Il989l : iHartmann & Kenyon"l996'; 'Reipurth & Aspin"2004'; 'Lorenze tti et al.ll2007l) . 

Significant advances in observational time domain astronomy are uncovering new 
phenomena and systems for study. This expansion is due to many factors: a steady im- 
provement in the size and sensitivity of optical arrays, enhancing the coverage and ca- 
dence possible for observations of optically revealed clusters; even greater advances in 
the capabilities of near-infrared arrays; and access to precise multi-epoch mid-infrared 
(mid-IR) photometry and spectroscopy from the Spitzer Space Telescope. These new 
observational capabilities have allowed variability studies to cover a greater number of 
targets at higher cadences, and characterize the variability properties of more deeply 
embedded, and presumably less evolutionarily advanced, sources. The signature of 
these advances is evident in our incr eased sensitivity to variabil ity at the lowest masses 
within nearby star-forming regions (ICodv & Hillenbrandll2010h as well as the increas- 
ing frequency with whi ch we id entify formerly rare young variables, incl uding new 
eclips ing binaries (.Hebb et al .*2010). new disk occulting systems (e.g., WL4: |Plavchan et al 
I2OO8I ). pulsating protosta rs (Morales-Calderon et al. 2009 ) and l arge-amphtude out - 
bursts dCovev et all I2OIOI : iMiller et all 12010 : iKospal et alTl2010l : iGaratti et all l2010h . 
These observational advances have been accompanied by similar progress on the theo- 
retical front, with increasingly detailed models of the physical proce sses underlying the 
observed variations (e.g., accretion variations, d isk processes, etc.; IVorobyov & Basu 
2OIOI : IZhu et al.ll2010l : iBaraffe & Chabrieill2010h 



'or, perhaps more relevantly, the timescale of a PhD thesis 
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This splinter session was convened to review recent progress in characterizing, an- 
alyzing, and understanding variability in the youngest stars. The session incorporated 
presentations covering physical processes that develop over a range of time-scales, with 
observational signatures spanning a wide range of wavelengths. Following these pre- 
sentations, the audience participated in a broader discussion of these new results and 
highlighted areas of particular promise for future observational or theoretical work. We 
provide a brief summary of each presentation below, and conclude with a recap of the 
open questions highlighted in the audience discussion. 



2. Mid-infrared (Mid-IR) variability 

2.1. Maria Morales-Calderon: A global view of Mid-IR variability from the 
YSOVAR Orion survey 

The YSOVAR Orion program provides the first large-scale survey of the photometric 
variability of Young Stellar Objects (YSOs) in the mid-IR. In Fall 2009, Spitzer/IRAC 
observed a 0.9 sq. deg. area centered on the Trapezium cluster twice a day for 40 con- 
secutive days, producing 3.6 and 4.5 jum high fidelity light curves (~3% typical pho- 
tometric uncertainties) for over 2000 disked and diskless YSOs in Orion (as diagnosed 
with precursor 3.6-8.0 yum IRAC photometry). For brevity, we refer to the diskless stars 
as Weak T Tauri stars (WTTs); we refer to stars with disks, or which are otherwise 
heavily embedded in their natal cloud, as YSOs. For many of the stars, we also ob- 
tained complementary time-series photometry at optical (1^) and/or near-infrared (JK,) 
wavelengths. We find that 65% of the disked YSOs and 30% of the diskless WTTs are 
variable. The WTT variations are mostly spot-like, while the disked/embedded YSOs 
display a wider range of variability types: see Figure [T] for examples of the different 
types of variability evident in the YSOVAR light curves. Consistent with their greater 
propensity for spot-like variability, 65% of the variable WTTs are periodic, while pe- 
riods are detected for only 16% of the YSOs. Those disked/embedded YSOs that are 
detected to be periodic, however, typically evidence larger amplitudes, as well as longer 
periods (see Figure |2ll. 

YSOVAR monitoring has also identified several rare pre-main sequence variables. 
Among these are five new candidate PMS eclipsing binaries, including one that is 
fainter than all previously known ONC PMS eclipsing binaries. One of the most sur- 
prising and interesting classes of variables we find are chara cterized by short du ration 
flux dips (less than one to a few days; ie, AA Tau analogs: iBouvier et al.ll2003h . We 



interpret these events as the star being extincted by either clouds of relatively higher 
opacity in the disk atmosphere, or geometric disk warps of relatively higher latitude 
which pass through the line of sight to the star. This hypothesis is consistent with both 
the duration and the wavelength dependence of the events (i.e., larger amplitudes at 
bluer wavelengths). 

2.2. James MuzeroUe: Spitzer Mid-Infrared Spectroscopic Variability in IC348 

We initiated a multi-epoch study of the young cluster IC 348 with all instruments on 
board Spitzer. This region is an appealing target given its reasonably close distance 
and wide range of YSO evolutionary states, including a significant number of highly 
evolved and transitional disks. With prior GTO/legacy observations, our dataset in- 
cludes a total of 10 maps of the entire cluster with MIPS, 7 maps with IRAC, and 4 to 7 
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Figure 1 . A sample of the range of variability types captured by the YSOVAR 
program. Most variations are relatively colorless (e.g., the bottom four panels), but 
some strong color variations are observed. Both red-dominant and blue-dominant 
color variations are detected (see top left and right panels, respectively), suggesting 
that different physical mechanisms (extinction, variable disk temperature & geom- 
etry, etc.) may be required to explain the full range of color-dependent variability. 
Aperiodic variability (which we typically interpret as disk/accretion driven) is the 
dominant mode of variability seen in the disked/embedded sources (e.g., middle-left 
panel), but we do see examples of periodic variability even in the most heavily em- 
bedded sources (e.g., middle-right panel). Unphased and phased light curves are 
shown for a particularly short period WTTs system (~0.27 d) in the bottom right 
and left panels, respectively. While we typically interpret periodic WTTs as exhibit- 
ing spot-induced variability, this system's short period and asymmetric light curve 
suggest it may be a close, tidally distorted binary. 
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Figure 2. 3.6 //m amplitude vs. period for YSOs detected as periodic variables 
in YSOVAR's mid-IR monitoring of Orion. Both period and amplitude appear cor- 
related with SED class; while periodic variability is most commonly detected for 
WTTs/Class III objects, those disked/embedded Class I-II YSOs that do demonstrate 
periodic variability have generally longer periods and larger amplitudes. 



epochs of IR S spectroscopy for 14 individual cluster members (InfraRed Spectrograph; 
Houck et"aD l2004). The measurements sample a range of cadences from days to years 
over the period 2004 - 2009. We find ubiquitous variability at all wavelengths: over 
half of the Class O/I objects at 3.6 - 24 fim, nearly 70% of Class II objects at 3.6 - 8 jum, 
and 40% of Class II objects at 24 /urn. Variability time-scales range from days to years 
at all wavelengths. We see apparent wavelength dependences with two general flavors: 
correlated variability over the 3.6 - 24 fim range, with or without color changes (most 
common in Class O/I sources), and anti-correlated behavior with a pivot point typically 
around 6-10 microns (most common in Class II sources and transitional disks). From 
the IRS spectra, almost all of the variations are related to the continuum; we see mea- 
surable changes in the silicate features in only a few objects. 

In most cases, the variations seen at longer wavelengths occur on time-scales too 
short to be representative of in-situ structural changes. For the Class II objects with 
anti-correlated wavelength behavior, the flux changes are likely connected to the in- 
nermost disk regions (puffed inner rim associated with dust sublimation, or dynamical 
truncation from close companions). Perturbations can cause a change in the inner disk 
height, which leads to direct changes in the short wavelength emission, and can cause 
shadowing of cooler material farther out in the disk, which leads to changes in the 
longer wavelength emission. The origin of the inner disk perturbations is not clear, 
but may be a result of changes in the accretion luminosity emitted by stellar accretion 
shocks, or inner disk warping created by tilted stellar magnetic fields or gravitational 
interactions with embedded companions. 

2.3. Kevin Flaherty: LRLL 31, a case study of Mid-IR variability in transition 
disks 

In addition to providing a robust characterization of mid-IR photometric variability, the 
Spitzer Space Telescope has diagnosed mid-IR variability in specific young stars via 
multi-epoch spectroscopy with IRS. These observations identified the detailed wave- 



6 



Covey et al. 



length dependence of mid-IR flux variations, pointing the way to new and unexpected 
results. One of the most interesting cases is that of LRLL 31, a G6 T Tauri star in 
IC 348: IRS spectra separated by one week revealed that the S-Sjum flux decreased to 



nearly photospheric levels while the 8-40//m flux increased by 60% (iMuzeroUe et al. 



20091) . What makes this star especially interesting is that it shows a deficit of flux 



around 10 //m compared to a traditional T Tauri star, indicative of a removal of small 
dust grains from the inner disk. This transition disk represents an evolutionary stage 
where the dust and gas are being removed from the circumstellar disk, and the fluctu- 
ations may be related to the process that is clearing the disk. Previous theories for ex- 
plaining variability in young stellar objects (rotation of hot/cold spots across the stellar 
surface, changes in the accretion rate, extinction events) do not explain the wavelength 
dependence or the strength of the variability observed in LLRL 31. Fluctuations on 
daily to weekly time-scales are surprising, since thermal equilibrium arguments would 
indicate that mid-IR wavelengths probe regions of the disk from 0.5 to a few AU, where 
the dynamical timescale is closer to years. In situ changes in this region of the disk 
should not be possible on these time-scales, so LLRL 3rs mid-IR variability must be 
rooted in the dynamics of the inner disk. 

Motivated by these observations of LLRL 31. lFlaherty & MuzeroUd (l20I0l) devel- 
oped a model that explains the observed variability as arising from a warp in LLRL 3 1 's 
disk, as opposed to the purely axisymmetric disks that had been previously assumed. 
Varying the height of a warp at the inner edge of the disk was able to explain the ob- 
served strength and wavelength dependence of the variability. As the warp grows it is 
more directly illuminated by the central star and heats up, emitting more flux at shorter 
wavelengths. The larger warp also shadows more of the outer disk, reducing the long- 
wavelength flux that originates far from the star. This model assumes the variability 
arises from structural changes in the disk warp: the precession of such a structure is 
unable to reproduce the variability. 

While this model was successful at explaining much of LLRL 3rs variability at 
A < 10/im, it was still limited by the fact that the 5-40 fim IRS spectra only partially 
traces emission from the inner disk. Observations at shorter wavelengths, produced by 
emission from the hottest dust, are needed to better understand the properties of the 
material responsible for the mid-infrared variations. Further monitoring at 3.6 and 4.5 
//m with Spitzer, and ground-based 0.8-5 jum spectra have revealed that the temperature 
of LLRL 3rs inner disk remained constant while the emitting area changed by a fac- 
tor of 10. Moreover, in 2009 there appeared to be a correlation between the accretion 
rate, as measured by the Pa^ line, and the infrared flux, measured by the infrared pho- 
tometry (Flaherty et al. submitted). This strong correlation suggests that the process 
that disturbs the dust must also disrupt the accretion flow. The most promising ex- 
planations are either a companion beyond 0.3 AU on an orbit misaligned with the disk 
(e.g. Fragner & Ne lson 200 9J or a dynamic interface between the s tellar magnetic field 
and the disk (e.g., ICoodson & Wingleelll999l : iBouvier et al.ll2007h . Contemporaneous 
X-ray and Spitzer observations will be able to trace some of the interaction between 
the magnetic field and the disk, further probing this process. The wavelength depen- 
dence of the variability seen in LRLL 3 1 appears to be common among transition disks 
(Espaillat et al. submitted) and further observations of both the gas and dust in these 
objects could determine if the physical process occurring in LRLL 31 is common in 
these other objects. 



Young Stars in the Time Domain Splinter Summary 



7 



2.4. Neal Turner: Models of axisymmetric disk disturbances from thermal waves 
and turbulent eddies 

T Tauri stars' diverse mid-IR variability likely arises from several different processes. 
Variations that correlate poorly with the accretion luminosity are most simply explained 
by changes in the way the starlight falls across the disk surface. Kevin Flaherty presents 
above some ways to produce non-axisymmetric disk disturbances. Among processes 
leading to approximately-axisymmetric disturbances, two explored here are thermal 
waves and turbulent eddies. 

Thermal waves are driven by the stellar illumination. The front of each wave tilts 
into the starlight, receiving extra heating and expanding, while the back of each wave 
receives less heating and cools and contracts. The wave propagates inward, dying out 
when it enters the shadow of the disk's inner rim. The wavelength is about equal to the 
length of the shadow cast by each peak. Thermal waves operate at mass accretion rates 
below about 10~' Mo/yr. For systems with larger accretion rates, accretion heating 
within the disk overwhelms the heating from external starlight, so that the waves no 
longer form or propagate. The wave pe riods are compa rable to the heating and cooling 
time-scales, and are measured in years (IWatanabe & L in 2008). 

But some T Tauri stars vary on time-scales of weeks, much shorter than the thermal 
wave period. The changes must arise near the disk inner edge, where the orbital period 
is about a week. On such short time-scales, the outer disk can be treated as static. A 
process producing changes in the disk photosphere height over times comparable to 
the orbital period is the same magneto-rotational turbulence that drives the accretion 
flow onto the star. The strength of the magnetic fields varies episodically, lead ing to 



vertical excursions of the disk photosphere by as much as 30% ( Turner et al.ll2010h . The 
excursions are approximately axisymmetric, at least over the length scales modeled in 
shearing-box calculations. When higher, the rim intercepts more starlight. Monte Carlo 
radiative transfer calculations (Turner N. J. et al., in preparation) show raising the rim 
increases the area with temperatures around 1000 K, while shadowing the disk beyond 
so that the emission from cooler material is reduced. The resulting variations in the two 
warm-Spitzer bands are correlated, with amplitudes generally decreasing toward longer 
wavelengths (see Figure (3]). Similar effects can be expected from the shadow cast by 
the outer edge of the gap opened in the disk by a giant planet, but the weaker frontal 
illumination would tend to decrease the amplitude. 



3. Short wavelength Variability 

3.1. Fabienne Bastien: High frequency variability of V1647 Orionis 

Numerous young stars have been observed to undergo large amplitude photometric out- 
bursts. While the observational record is still woefully incomplete, these events may 
play a significant role in the early evolution of many, possibly even most, young stars. 
These outbursts have traditionally been divided into two major classes: FU Ori-like 
large amplitude (> 4 mag), long duration (t> 10s years) outbursts, w ith spectra domi- 



nated by absorption features characteristic of cool (GKM) supergiants (iHartmann & Kenyon 
Il99 6l). and EX Lup-like (often shortened to EXor) outbursts, which appear to be more 
frequent, of shorter duration, and characterized by heavily veiled spectra with strong 
emiss ion and absorption components in accretion and outflow-sensitive features dHerbigl 
l2007h . 
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Figure 3. (Left) - A bump near the inner rim casts a long, cool shadow across the 
disk while itself being heated by the intercepted starlight. Curves show the surfaces 
of unit optical depth for light coming from the star (solid) and from directly overhead 
(dashed), for cases with (red) and without a bump (green). (Right) - Bumps whose 
height changes over time produce variability with amplitude greatest at shorter in- 
frared wavelengths. Synthetic IRAC 3.6 and 4.5 jum light curves for a model disk 
with a magnetically inflated inner rim. (both panels. Turner et al., in prep). 



In 2004, the star VI 647 Orionis was observed to undergo a major accretion out- 
burst. Numerous observations were obtained throughout the outburst, which persisted 
through 2006, but the source's photometric and spectroscopic properties resisted simp le 
classi fication as a prototypical FU Ori-like or EX Lup-like outburst (e.g., iFedele et all 



20071 ). Archival photometry revealed that the 2003 event was not V1647 Ori's first 
outburst, and additional obje cts have been obse rved to undergo outbursts with similar 
spectral characteristics (e.g.. ICovey et al.ll2010l) . suggesting that these outbursts may 



signify a new phenomena of import for a star's early evolution. 

Our study of high cadence time-series photometry of the 2003-2004 and 2008- 
2009 eruptions of VI 647 Orionis revealed interesting transient variability on time- 
scales of a few hours: we detect a highly significant period of 0.13 days in our 2003 
data that does not appear in our 2009 data. We attribute this period, which is inconsis- 
tent with the object's rotation period, to a short-term radial pulsational mode of the star 
excited by the sudden increase in its accretion rate at the time of our 2003 observations. 
This period, if also excited in VI 647 Ori's latest outburst event, would have presumably 
rung down by the time we observed the star in 2009. More near infrared and optical 
observations with cadences of at least one image per hour over the course of several 
days during the rise, plateau and decline phases of such outbursts would allow us to 
better understand what short timescale phenomena occur during these events. 

3.2. Ettore Flaccomio: Correlated optical and X-ray variability in accreting 
Classical T Tauri Stars 

The magnetospheric accretion model, in which material is funneled from the inner 
edge of the circumstellar d isk to the stellar photosphere along magnetic field lines 



(e.g., lOstriker & Shulll995l : iBouvier et al.ll2007l) . predicts that accreting material will 



be strongly shoc ked when it impacts the stellar surface. This extremely hot (T ~ 10,000 



K: lHerczeg & H illenbrand 2008) shocked region is believed to pr oduce significant con 



tinuum emission at optical, ultraviolet, and X-ray wavelengths dGiinther et al.ll20()7 ). 



Clear differences have indeed been detected between the X-ray properties of accret- 
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ing Classic T Tauri stars (CTTs) and non -accreting WTTs: CTT s possess harder, more 
time-variable X-ray spectra than WTTs dFlaccomio et al.ll2006l) . and detailed spectral 
studies sugges t that accretion con tributes a distinct, soft (2-3 MK) component to CTTs 
X-ray spectra dOudel et al.ll2007h . 

Many open questions remain, however, as to the exact mechanisms underlying the 
connection between accretion and X-ray emission. For instance, accreting systems are 
counter-intuitively seen to be significan tly less X-ray lum inous than predicted by stan- 
dard magnetospheric accretion models dJohns-Krul l' 2007) as well as their WTTs coun- 
terparts (P reibisch et al..2005) . potentially due t o the shocked emissio n being extincted 
by material in the overlying accretion columns ( Gregory et al.ll20()7h . Previous obser- 
vations have also failed to identify a clear correlation between CTTs optical and X-ray 
emission. Accretion is expected to contribute a portion of the X-ray e mission, origi- 
nating in the photospheric accretion shock (as opposed to coronal X-rays: [Stassun et aL 
200al2007l:rGrosso et al.ll2007h. 



In March 2008, lAlencar et all fcOlOh observed NGC 2264 with CoRoT for 23.5 
days obtaining high-quality uninterrupted optical light-curves of its young stars. During 
the CoRoT pointing, two short Chandra observations were performed with a separation 
of 16 days, allowing us to study the correlation between optical and X-ray variabil- 
ity on this timescale, a nd thus the physical mechanism responsible for the variability 
dFlaccomio et al.ll2010l) . The variabilities of CTTs in the optical and soft X-ray (0.5-1.5 
keV) bands are correlated, while no correlation is apparent in the hard (1.5-8.0 keV) 
band. Also, no correlation in either band is present for WTTs. The correlation be- 
tween soft X-ray and optical variability of CTTs can be naturally explained in terms 
of time-variable shading (absorption) from circumstellar material orbiting the star, in a 
scenario rather similar to the one invoked to explain the observed phenomenology in 
the CTT star AA Tau. The slope of the observed correlation implies (in the hypothesis 
of homogeneous shading) a significant dust depletion in the circumstellar material. 



3.3. Marsden 

One of the key aspects of the solar dynamo is the reversal of the Sun's global magnetic 
field every ~11 years. However, the dynamo mechanism for young rapidly -rotating 
solar-type stars is unclear and may well be fundamentally different to that of today's 
Sun (Donati et al. 2003 ). To learn more about the type of dynamo opera ting in young 
stars we have used th e technique of Zeeman Doppler Imaging (ZDI; ISemel|]l989l : 
iDonati & Browniri997l) to map the magnetic topology of two young solar-type stars 



(HD 141943 and HR 1817) over several years, searching in particular for evidence of 
polarity revers als in the global mag netic field. Given the rapid differential rotation seen 
on such stars (Marsden et alllOlOah it is expected that they have significantly shorter 
magnetic cycles than the Sun, perhaps as short as a few years. 

While HD 141943 shows evidence of changes in its magneti c topology it shows n o 



evidence of a polarity reversal over the 3 years of observations dMarsden et aLllOlOah . 



This result is similar to another young solar- t ype star (HD 171488 ) which has been 
previ o usly observed usin g similar techniques dMarsden et all l2006l : IJeffers & Donati 



20081 : Ijeffers et al.l 120101) . In contrast, the preliminary results for HR 1817 indicate 



a reversal of the polarity of the star's radial magnetic field between 2008 and 2009, 
although a similar reversal is not seen in its azimuthal magnetic field (see Figure |4l 
Marsden et al., in prep.). Such behavior suggests a phase delay in polarity reversal 
between the radial and azimuthal fields. 
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Figure 4. Preliminary radial (top) and azimuthal (bottom) magnetic field topolo- 
gies for HR 1817 for the years (left to right), 2001, 2007, 2008 and 2009. The images 
are flattened polar projections extending down to -30° latitude. The bold line denotes 
the equator and the dashed lines are -1-30° and -1-60° latitude parallels. The radial ticks 
outside the plot indicate the phases at which the star was observed and the scale is in 
Gauss. (Marsden et al., in prep.) 

In order to assess how young solar-type stars undergo magnetic polarity reversals, 
optical spectropolarimetric observations of a statistically useful number of such stars 
needs to be undertaken regularly over a number years. A determination of how young 
solar-type stars go through magnetic polarity reversals and whether they have regular or 
chaotic magnetic cycles should greatly help our understanding of the dynamo in young 
Suns. 

4. Conclusions 

As this splinter summary demonstrates, recent years have seen considerable progress in 
our ability to observe, characterize, and model time- variable processes in the formation 
and early evolution of stars, circumstellar disks and planets. These advances include 
the ability to study variability across a wider range of wavelengths and time-scales 
than previously possible, and to construct more detailed and computationally intensive 
models of star and disk processes. 

Nonetheless, our current understanding of time-variable phenomena in early stel- 
lar evolution is significantly incomplete. In addition to the questions highlighted above, 
the audience identified several other outstanding challenges for studies of time-domain 
processes in the formation and early evolution of stars and planets. These include: 

• Improving the mechanisms and capabilities for conducting multi-site and multi- 
wavelength studies with long time baselines. 

• Achieving statistically robust constraints on variability (e.g., FU Ori outburst 
rate/duty cycle), rather than qualitative descriptions. 

• Identifying a critical mass of individual, rare systems (i.e, KH 15D, McNeil's 
Nebula, etc.) such that we can infer universal lessons from their specific proper- 
ties. 
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• Developing better observational and theoretical diagnostics to locate the source 
regions and determine the causes of variability. Two examples are detecting kine- 
matic signatures of the underlying gas motions in the line emission from the star 
and disk, and linking the photometric variability of the inner disk to changes in 
the scattered light from the spatially-resolved outer disk. 

• Understanding the implications for planet formation of time- variable structure in 
protoplanetary disks. 

Meeting each of these challenges will require ingenuity, focused effort, and careful 
planning. The young stars presenting here today suggest that we will indeed be up to 
the task, but only tim^l will tell. 

Acknowledgments. We thank the Cool Stars 16 SOC for organizing a productive 
and enjoyable meeting. We particularly thank the editors of this volume, whose enthu- 
siasm and encouragement were critical to the "timely" preparation of this article. 
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